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9 Abstract 
10 

1 1 When averaged over the tropical oceans (30° N/S), latent heat flux anomalies derived from passive 

12 microwave satellite measurements as well as reanalyses and climate models driven with specified sea- 

13 surface temperatures show considerable disagreement in their decadal trends. These estimates range 

14 from virtually no trend to values over 8.4 Wm" 2 decade' 1 . Satellite estimates also tend to have a larger 

15 interannual signal related to El Nino / Southern Oscillation (ENSO) events than do reanalyses or model 

16 simulations. An analysis of wind speed and humidity going into bulk aerodynamic calculations used to 

17 derive these fluxes reveals several error sources. Among these are apparent remaining intercalibration 

1 8 issues affecting passive microwave satellite 10 m wind speeds and systematic biases in retrieval of 

1 9 near-surface humidity. Likewise, reanalyses suffer from discontinuities in availability of assimilated 

20 data that affect near surface meteorological variables. The results strongly suggest that current latent 

21 heat flux trends are overestimated. 
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Introduction 


2 

3 Latent heat flux (LHF) over the global oceans is a key climate process linking fast convective, 

4 radiative, and dynamical processes of the atmosphere to the slower, but equally dynamically active 

5 thermal reservoir below. Prospects for enhancing climate predictability from seasonal to decadal scale 

6 and beyond depend critically on an improved understanding and modeling of physical processes 

7 controlling LHF. The WMO Working Group on Air-sea Fluxes [Taylor, 2000] has articulated this 

8 need in broadest terms. The SEAFLUX project [Curry et al., 2004] has emerged as an example of 

9 focused efforts to evaluate and improve satellite-based products in particular. As a result of these 

10 efforts estimates of LHF variability on scales from interannual to decadal are now emerging. In 

1 1 particular, evidence for upward trends in LHF over the tropical domain has recently been reported. 

12 Chiu and Xing [2004] performed an EOF analysis on version 2 of the Goddard Satellite-based Surface 

13 Turbulent Fluxes data set (GSSTF2) and found that the leading mode exhibited an upward trend over 

14 the period 1990-2000. The structure of this mode showed positive changes in subtropical dry regions 

15 and smaller negative changes in the equatorial Pacific. Liu and Curry [2006] expanded this analysis to 

16 include reanalyses and an additional satellite algorithm, HOAPS-II (Hamburg Ocean Atmosphere 

17 Parameters and Fluxes from Space version II). Their analysis provided a similar finding of a 

1 8 statistically significant upward LHF trend having a similar spatial structure and they argued that 

1 9 surface wind speed increases were a key driver. It should be noted that the rate of LHF increase 

20 inferred from the GSSTF2 by both these studies was about 17 Wm' dec' , while the reanalysis 

21 estimates in Liu and Curry range from 1.6 to 8.9 Wm' 2 dec'\ Given the estimated climatological 

22 tropical mean LHF value of ~ 100 Wm' 2 , variations of this magnitude demand further scrutiny. In this 

23 paper we re-examine several data sets and attempt to understand and narrow the uncertainties. 

24 
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Data and Methodology 


1 2 . 

2 

3 The analysis presented here relies on monthly means of LHF estimated through bulk aerodynamic 

4 methods as applied to satellite data and in numerical models. Two satellite-derived products are used, 

5 both of which rely heavily on Special Sensor Microwave Imager (SSM/I) radiances to provide 10m 

6 wind speeds and estimates of near surface specific humidity, q a . For both of these data sets the fluxes 

7 were derived on daily time scales and archived as monthly mean products. HOAPS-II provides ocean 

8 turbulent flux retrievals and related quantities on several spatial and temporal scales [Schultz et al., 

9 1997; Grassl et al., 2000; Bentamy et ah, 2003]. Sea-surface temperatures (SSTs) taken from the 

1 0 Advanced Very High Resolution Radiometer Pathfinder data set are used to compute saturation 

1 1 specific humidity at the ocean surface, q s . We used the monthly mean, 0.5 degree lat / Ion resolution 

12 data for the period Jan. 1988 through December 2002, omitting the period July 1991 through July 1992 

13 when aerosols from the Pinatubo eruption contaminated SST retrievals. A second data set from the 

14 University of California, Santa Barbara (UCSB) providing 0.25 degree resolution monthly mean 

15 retrievals was also used. The UCSB retrieval [Jones et al., 1999] uses a neural net methodology, 

16 training on the NOAA/PMEL TAO and Pirata buoy arrays [McPhaden, 1995] to retrieve 2m q a (and 

17 T a ). Wind speeds for this algorithm are the version 5 retrievals of Wentz (1997). SSTs are taken from 

1 8 the NCEP/NCAR (National Center for Environmental Prediction / National Center for Atmospheric 

19 Research) Reanalysis [Kalnay et al., 1996]. 

20 

21 Monthly mean fluxes and near surface meteorology were also taken from three reanalyses. These 

22 included the NCEP/NCAR Reanalysis [Kalnay et al., 1996], the NCEP-2 Reanalysis [Kanamitsu et al., 

23 2002], and the new Japanese Meteorological Agency Reanalysis, JRA-25, [Onogi et al. 2005]. These 

24 data sets will be referred to as NN, N2, and JRA-25, respectively. The JRA-25 assimilated SSM/I 

25 column-integrated water vapor (CWV) as well as ERS-1, 2 and QuikScat scatterometer wind speeds. 
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1 However, neither the NN nor the N2 use either of these, thus providing a useful measure of 

2 independence from the other data sets. Finally, we include LHF produced by an integration of the 

3 NASA FVGCM model [Lin, 2004] forced by observed SSTs. The data are from the same integration 

4 analyzed by Bosilovich et al., [2005]. This integration is also a useful comparison because even 

5 though it is constrained by the time variation of SSTs, it is free from any time-dependent biases 

6 associated with data assimilation (e.g. the availability of SSM/I data only after July 1987). The 

7 physical parameterizations of the FVGCM are based on National Center for Atmospheric Research 

8 (NCAR) Community Climate Model version 3.0 (CCM3) physics. 

9 

10 3. Analysis 

11 

12 Our analysis focuses on tropical ocean area-averages, more specifically on departures from 

13 climatological monthly means over the base period January 1988 to December 2002. Figure 1 contains 

14 time series of LHF, q s *- q a , and 10 m wind speed anomalies area-averaged over the tropical oceans 

15 (30° N/S) for the data sets described above. Also shown in shading are associated SST anomalies 

16 calculated in the same fashion. Each time series has been smoothed using a 5-month running mean for 

17 presentation clarity. The LHF time series (Fig la) each show increases with time, although there is 

18 substantial variability, ranging from 2.6 Wm" 2 decade _1 for the FVGCM to 8.4 Wm‘ 2 decade _1 for UCSB. 

19 This unanimity in trend might seem impressive were it not for the fact that these values are several 

20 time larger than GPCP tropical oceanic rainfall trends (1.2 Wm' 2 decade _I ). Since atmospheric water 

21 vapor storage is vanishingly small on these scales, there would be a large implied atmospheric 

22 transport of excess moisture to support increased precipitation over tropical land or higher latitudes. 

23 This is not apparent in the GPCP record which actually shows small rainfall decreases over land areas 

24 with time [Gu et al., 2007]. Note also that the FVGCM, though capturing much of the interannual 

25 variability associated with ENSO, has the smallest trend. While the accuracy of atmospheric model 
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1 integrations forced by observed SSTs is open to question [Wang et al., 2005], one must also question 

2 the effects of data set inhomogeneity in assimilated fields, and the accuracy of satellite-based LHF 

3 algorithms. 

4 

5 The origin of these differences becomes somewhat clearer when we examine time series of q s *- q a (Fig 

6 lb). Of all the model runs, JRA25 shows the greatest rate of increase (0.21 g Kg" 1 decade' 1 ). Part of 

7 this can be attributed to the effects of the onset of SSM/I data availability in July 1987. Relative to NN 

8 and N2 which did not assimilate SSM/I data, the JRA25 is lower by order 0.10 g Kg' 1 before this date 

9 as compared to the four years following. It is also apparent that N2 anomalies exhibit a larger trend 

10 than NN while having similar interannual variability. Before 1995 the offset N2 anomalies are larger 

1 1 negative, but that after this date N2 anomalies increase to values nearly 0. 10 g Kg' 1 larger than NN. 

12 SST and, thus, q s * is essentially invariant between these two reanalyses. In contrast, the FVGCM 

13 integration has a much smaller q s - q a trend (0.04 Wm' decade' ). 

14 

1 5 Another interesting aspect in Fig. lb is the distinctly larger magnitude of q s *- q a anomalies in the two 

1 6 satellite-based algorithms during ENSO events compared to those in the reanalyses. We believe that 

17 this stems from the inability of satellite algorithms to discern systematic variations in the relationship 

18 of near-surface moisture in the planetary boundary layer (PBL) with that of CWV. This issue was first 

19 noted in a climatological sense by Esbensen et al. [1993] who found that Liu’s (1986) method of 

20 retrieving q a from CWV yielded systematic underestimates of q a in subsiding trade wind regimes and 

21 over-estimates in zones of convergence and deep convection. In the subtropics, mean tropospheric 

22 subsidence in the atmospheric column dries the column by vertical advection. However, in the well- 

23 mixed lower portion of the PBL, the fractional reduction in moisture is not as pronounced as in the 

24 free troposphere. In regions of anomalous ascent, vertical pumping of moisture provides a larger 

25 fractional elevation of free-tropospheric moisture than occurs near the surface. To see if this effect is 
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1 modulated by ENSO we constructed satellite minus NN q a anomalies and composited them for both El 

2 Nino and La Nina years. By further differencing the La Nina from the El Nino composites (e.g. warm 

3 SST versus cold SST years) we see systematic patterns of q a bias in UCSB or HOAPS q a relative to 

4 NN (Fig. 2a, b). Overestimates in the central equatorial Pacific during El Nino are surrounded by 

5 underestimates to the west and over the adjacent subtropics. This pattern corresponds to anomalous 

6 upward motion in the equatorial eastern Pacific and sinking motion in the surrounding areas. Bias 

7 values are more prominent in HOAPS than in UCSB. The effects of widespread tropospheric sinking 

8 during warm events dominates the area average value in both data sets, leading to overestimates of q s *- 

9 q a and, thus, LHF anomalies. The opposite situation holds during ENSO cold events (e.g. 1999-2000). 
10 

1 1 Another way of looking at this issue is to compare the relative variations of CWV and q a from 

12 reanalyses during ENSO events. All q a estimates based on satellite passive microwave rely essentially 

13 on the CWV signal either explicitly or implicitly. Thus, for consistency with the bias mechanisms 

14 outlined above we would expect that true fractional CWV variations are larger than those of true q a . 

15 Using NN data we composited anomalies of logio CWV and logio q a for El Nino and La Nina events, 

16 subtracting the La Nina composite from the El Nino composite. The results (Fig 2c, d) indeed show 

17 that compared to q a , fractional increases in CWV are larger positive over the equatorial Pacific and 

1 8 larger negative over subsidence regions in the subtropics and Maritime Continent. This is not a result 

1 9 peculiar to NN, the same results were obtained when we used the FVGCM which is independent of 

20 any data assimilation uncertainties. On the basis of these results we conclude that the pronounced 

21 variations in q s *- q a for the satellite data sets (Fig lb) are likely exaggerated, as are interannual 

22 variations in LHF (Fig la) that they support. 

23 

24 Figure lc compares monthly wind speed anomalies (computed from daily wind speeds). The 

25 magnitudes of these departures are several tenths of m s' 1 corresponding to order 5% departures from 
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1 climatological values. NN and N2 anomalies track very closely, perhaps because assimilated 

2 dynamical fields are more strongly shielded from model physics changes than are quantities such as 

3 near-surface moisture. Likewise, JRA-25 tracks NN and N2 before 1995, but then decreases for 

4 several years before suddenly rising sharply in 2001. Mid- 1995 marks the onset of Environmental 

5 Research Satellite (ERS-1, 2) scatterometer wind assimilation by JRA-25. QuikScat winds are 

6 assimilated beginning in mid- 1999 and are the only scatterometer winds after ERS-2 ceases in 2001 . It 

7 is possible that rain and cloud contaminated wind retrievals may still not be adequately flagged in these 

8 data sets and are having a systematic effect. Interestingly, the FVGCM wind speeds have little 

9 coherence with reanalysis winds on short time scales associated with transient climate “noise” but 

10 seem to capture lower frequency behavior such as elevated speeds in the late 1980s and lower values in 

1 1 the early 1990s. They also show no signs of the JRA-25 wind speed drop and subsequent rise 

12 beginning in 1995. 

13 

14 Perhaps the most surprising aspect of Fig.lc is the large variability and trend associated with the two 

15 SSM/I wind time series. These algorithms are independent from each other yet agree closely in their 

16 departures from assimilated wind speeds. Here too it is possible that liquid water is contaminating 

17 wind retrievals, but satellite intercalibration seems a more likely explanation. The agreement of high 

1 8 frequency interannual SSM/I variations with assimilated wind speeds (e.g. the mid-1990s) and the 

19 presence of step-like departures at the end of 1992 and 1997 support this argument. Wentz et al., 

20 [2007] have recently shown that intercalibration improvements in the new version 6 of the SSM/I 

21 products greatly reduce this trend in SSM/I surface wind speeds and substantially reduce inferred LHF 

22 estimates. 

23 

24 4. Discussion 

25 
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1 From this comparison of LHF data sets on interannual to decadal scales there clearly exist uncertainties 

2 in current data sets that rival interannual to decadal signals. The most prominent distinction is between 

3 satellite-based and reanalysis LHF estimates. The former appear to have significant explainable biases 

4 in q s *-q a on ENSO time scales, and the wind speeds have significant departures from assimilated fields 

5 suggestive of sensor intercalibration issues. The latter drive large upward temporal trends in the 

6 satellite LHF. The GSSTF2 data set analyzed by Liu and Curry [2006] has an even larger trend (~1 7 

7 Wm 2 decade' 1 ) due in part to use of the earlier version 4 of SSM/I wind and CWV products. The 

8 reanalysis with the largest trend, JRA-25, owes part of this behavior to non-uniform temporal extent of 

9 assimilated data sets, most notably, onset of SSM/I availability in 1987 and the use of scatterometer 

10 winds beginning in 1995. These results do not, however, suggest that NN and N2 should be regarded 

11 as definitive since there are many kn own deficiencies in these and other reanalyses. We have 

12 purposely excluded the 40-year European Center Reanalysis, ERA-40, because of satellite radiance 

13 calibration and water vapor assimilation issues which have led to known biases in hydrologic cycle 

14 processes [Bengtsson et al., 2004; Andersson et al, 2005]. It appears also that blended products using 

15 satellite, reanalysis, and buoys cannot fully resolve the issues found here. Yu and Weller [2007] report 

16 an increase of about 9 Wm' 2 globally during the period 198 1 to 2002 for their Objectively Analyzed 

17 Air-sea Heat fluxes (OAFlux) data set. However, the GSSTF2 q a estimates and winds taken from 

18 version 4 of the Wentz SSM/I products are integral components of this analysis and likely contribute to 

1 9 a spurious upward trend. 

20 

21 We would argue that given the error sources identified in the present results, the degree to which 

22 interannual variations in LHF at regional scales integrate to a discemable tropic-wide trend remains an 

23 open question. In terms of long-term greenhouse gas forcing effects, Held and Soden [2006] argue for 

24 a scaling of order 2% rise in precipitation (or evaporation) per degree SST rise which, given the -.17 K 

25 SST rise over the tropical oceans during the past 25 years, would yield a LHF increase of order 0.5 
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1 


2 ■ * 

Wm’ . This is far smaller than variations noted in the present results. However, there do exist 

2 candidate physical mechanisms for producing variations of decadal scale and longer. Among these are 

3 possible reddening of the ocean upper ocean heat content variations by “higher frequency” ENSO 

4 events similar to that hypothesized by Newman et ah, [2003] or low-frequency variations in the 

5 shallow meridional overturning circulation of the Pacific Ocean [McPhaden and Zhang, 2004]. 

6 Associated ocean heat transports would likely lead to subsequent adjustments with the overlying 

7 atmosphere through surface fluxes. 

8 

9 It is expected that ongoing improvements to existing data sets [e.g. Wentz et al., 2007] and planned 

10 reanalyses (e.g. the NASA Modem-Era Retrospective Analysis for Research and Applications) will 

1 1 substantially improve prospects for isolating physical signals of LHF variations. 

12 
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1 Figure Captions 

2 

3 Figure 1. Anomaly time series for LHF, q s *-q a , and 10m wind speed for NN (black), N2 (black dot), 

4 JRA-25 (blue), FVGCM (gold), FIOAPS II (green) and UCSB (red). Anomalies are with respect to 

5 1988-2002 mean annual cycle (except FVGCM, 1988-1998). SST anomalies are shaded gray. 

6 FVGCM wind speeds are only available through 1998. 

7 

8 Figure 2. Anomalies of (a) q a (UCSB minus NN) and (b) q a (HOAPS II minus NN) for El Nino minus 

9 La Nina composites with units of g Kg' 1 ; (c) fractional anomalies (dimensionless) of logio CWV and 

10 (d) logio 9a for El Nino minus La Nina composites made from NN reanalysis data. 

11 
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